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Stellar populations in NGC 5128 with the VLT: evidence for 

recent star formation* 
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Abstract. We resolve stars of the nearest giant elliptical galaxy NGC 5128 using VLT with F0RS1 and ISAAC. 
We construct deep U, V and K s color-magnitude and color-color diagrams in two different halo fields (in the halo 
and in the north-eastern diffuse shell). In the outer, shell field, at ~ 14 kpc from the center of the galaxy, there 
is a significant recent star formation with stars as young as 10 Myr, approximately aligned with the prominent 
radio and x-ray jet from the nucleus of the host AGN. Ionized gas filaments are evident in ultraviolet images near 
the area where neutral Hi and CO molecular gas was previously observed. The underlying stellar population of 
the halo of the giant elliptical is predominantly old with a very broad metallicity distribution. The presence of 
an extended giant branch reaching Mbol = —5 mag suggests the existence of a significant intermediate-age AGB 
population in the halo of this galaxy. 
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y-< 1. Introduction 



Galaxies show a wide range of star formation activity and 
a large range of mctallicities and ages in their stellar pop- 
ulations. Understanding the physical nature and origins of 
their stars is fundamental to understanding the formation 
and evolution of galaxies. The current knowledge of star 
formation histories of galaxies along the Hubble sequence 
is summarized by Kennicutt (1998). This knowledge is 
mainly based on the integrated photometric and spectro- 
scopic studies through the predictions of populations syn- 
thesis of the integrated light of star clusters (Bica 1988) 
and galaxies (Bruzual & Chariot 1993[ Maraston 1998[ ). 
The predictions of population synthesis can now be tested 
directly not only for globular clusters (e.g. Vazdekis et 
al. 2001), but also for nearby galaxies. 



Stellar evolution theory provides predictions of the 
features expected in color-magnitude diagrams (CMDs) 
for stellar populations with different ages and metal- 
licities (Renz ini & Fu si Pe cci 1988, Chiosi et al 



1992 



Gallart [L99f, Aparicio |i~99S| , Tolstoy |l998|) . Coupled with 
the improvements in telescope sizes, detector sensitivity, 
field of view and spatial resolution, compared with those 
a decade ago, the direct observations of the stellar con- 
tent of nearby galaxies is becoming one of the most active 
areas of extragalactic research. 

In the Local Group, indicators of the old stellar popu- 
lations such as old main sequence turn-offs or at least the 
horizontal branch magnitudes are within reach of available 
(e.g. Hurley-Keller et 
200C| . 
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NAS5-26555. 



recent results on the studies of the Local Group galax- 
ies have been summarized by Mateo (199S), van 
Bergh dl999| , p000| ) and Grebel ( p000| ). 



den 



The Local Group contains galaxies representative of 
almost all the classes, except the important giant elliptical 
class of galaxies. The closest giant elliptical is NGC 5128, 
the dominant galaxy in the Centaurus group at dist ance of 

Soria et al. 



3.6 ±0.3 Mpc (Harris et al. [1999 
al. |1995|, Tonry & Schechter 



1996, Hui et 



199C ) . It has been extensively 
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studied over the last 50 years (for an exhaustive review see 
Israel |1998| ). 

Its popularity is not only due to its brightness, but 
also to its many unusual features: (i) there is a promi- 
nent dust band containing young stars and Hn regions 
(Unger et al. |200(i| , Wild & Eckart p000| , Graham |l97§ ); 
(ii) there is an active nucleus with a radio and x-ray 
jet, radio lobes and optical filaments (Cooper et al. 
Feigelson et al. 1981, Kraft et al. |2000 , Schreier et al 



1965 



Clarke et al. 1992, Blanco et al. 1975. Dufour & van den 



1981 



Bergh |1978| ); (iii) Malin et al. fl!983| j discovered a large 



number of faint narrow shells o f star s surrounding the 
galaxy; (iv) Schiminovich et al. ( |l994[) detected 4 x 10 8 
M Q of Hi gas associated with the stellar shells, but slightly 
displaced outside the shells; and (v) most recently, molec- 
ular CO gas has been found associated wit h the Hi gas 
and the stellar shells (Charmandaris et al. 2000 ). All of 
these are clear indications of a recent interaction with a 
gas-rich galaxy. 

The high resolution and sensitivity of WFPC2 on the 
Hubble Space Telescope (HST) enabled the first stud- 
ies of the resolved old stellar populations in the halo of 
NGC 5128 (Soria et al. |l996j Harr is et al. |1999| Harris 
& Harris |2000[ Mould et al. [2000D . NICMOS on HST 
was used to resolve the stars in the near IR (Marleau et 
al. |2000|) in the same field as the optical study of Soria 



et al. (1996; ^9kpc south of the center of the galaxy). In 
these two studies, a small intermediate- age population of 
~ 5Gyr has been found, comprising up to ~ 10% of the to- 
tal stellar population in the halo. On the other hand, there 
are almost no intermediate-age stars in the field further 
out in the h alo, at ~ 20 kpc from the galaxy ce nter ( Harris 
et al. |1999[) nor at - 31 kpc (Harris & Harris |200C| ). The 



comparison of the two results may indicate the presence 
of a gradient in the stellar popu lation in the halo, as sug- 
gested by Marleau et al. (200C). However, the small field 
of view of HST puts serious limitations to the conclusions 
in the cases where the strong gradients in galaxy popula- 
tions exist (see, for example, the ca se of Local Group dwar f 
galaxies like Leo I (Gallart e t al. |1999| Held et al. [2000D 
or WLM (Minniti & Zijlstra |l996Hl997| )). 

We present here a wide wavelength range photometry 
of the resolved stellar populations in NGC 5128 obtained 
from the ground with Very Large Telescope (VLT) in U, V 
and K s band. The deep and high-resolution VLT imaging, 
coupled with the much larger field of view than HST, en- 
ables us to address the questions of the gradients in stellar 
populations in the halo of this giant elliptical galaxy. The 
proximity of NGC 5128 provides an unusual opportunity 
for a direct study of shell stars. We use infrared-optical 
colour-magnitude diagrams of the shell to study the ages 
and metallicities of the stars belonging the cannibalized 
galaxy. 

2. The Data 



The observations were carried out with the ESO Very 
Large Telescope (VLT) at Paranal Observatory, Chile. 



They consist of optical (Bessel U - and V-band) and near- 
IR (K s -b&nd) images. We observed two fields in the halo 
of NGC 5128. Field 1 (a 20 oo = 13 h 26 m 23.5 s , <5 20 oo = 
-42°52'00"; Fig. 0) was centered on the prominent N- 
E shell, ~ 14' away from the center of the galaxy. Field 2 
(a 20 oo = 13' 1 25™26 S , S 20 oo = -43°10'00"; Fig. Jwas cho- 
sen to overlap with the HST field of Soria et al. ( |1996 ) and 
lies at a distance of ~ 9' from the center of the galaxy. 
The optical data were obtained on 1999 July 11 and 12, 
while the if s -band images cover several epochs and will 
be used to search for long-period variable stars. All of the 
observations were secured in the service mode. Calibration 
data, bias, dark and flat-field images as well as photomet- 
ric stand ards fr om Landolt (in optical; 1992| ) and Persson 
et al. (IR; 199?) catalogues, were supplied by the ESO cali- 
bration plan. The observations are summarized in Tabic [l]. 

2.1. Optical Observations and Data Reduction 

Two 15-minute exposures in both U and V were ac- 
quired for each of the two fields in NGC 5128 using 
FORS1 (FOcal Reducer and low dispersion Spectrograph) 
on VLT Unit Telescope 1 (Antu). The FORS1 detector is 
a 2048x2048 Tektronix CCD, thinned and anti-reflection 
coated. The pixel size is 24x24 /xm. The field of view is 
6.'8 x 6f8 and the scale is O'.'2/pixel. 

For service observations in direct-imaging mode, the 
FORS1 CCD is read out in four-port read-out mode. We 
used the ESO pipeline reductions which are based on 
MIDAS package specially developed to reduce images with 
4 different amplifiers. With it the overscan for each ampli- 
fier was subtracted individually and after the subtraction 
of bias the images were corrected for the flat-field. 

2.2. Near-IR Observations and Data Reduction 

For the near-IR K s -band observations we used ISAAC, 
also on Antu. In this wavelength domain (0.9-2.5 /im) the 
detector is a 1024x1024 Hawaii Rockwell array. The field 
of view is 2.'5 x 2'.5 and the scale is O'.'147/pixel. A series 
of 10-second exposures were taken at each epoch, usually 
in groups of six, with the number of repeats depending on 
weather conditions and any technical problems. The total 
exposure time for each epoch is given in Table |. 

An additional if s -band observation of Field 2, with 
a total on-target exposure time of 45 minutes, was ac- 
quired with the SOFI instrument at the ESO 3.5-m New 
Technology Telescope (NTT) at La Silla Observatory, 
Chile. The field of view of SOFI is 4f94 x 4.'94 and the 
scale is O'.'292/pixel. 

We also obtained observations with the NIC3 array on 
HST, using the F222M filter, which is similar to (although 
narrower than) the standard K filter. The NIC3 field of 
view is 51'.' 2 x 51'.' 2, much smaller than the one of ISAAC. 

The standard procedure in reducing IR data consists of 
(i) dark subtraction, (ii) flat-field correction, (iii) sky sub- 
traction, (iv) registering and combining the images. We 
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Field 


a (2000) 


^(2000) 


Date 


Telescope & 


Exposure 


FILTER 


Airmass 


Seeing 


Epoch 


# 


(h 


min sec) 


(° ' ") 


dd/mm/yyyy 


Instrument 


(sec) 










1 


13 


26 


23.5 


-42 52 00 


12/07/1999 


Antu- 


-FORS1 


2x900 


U 


1.54 


0.52 




1 


13 


26 


23.5 


-42 52 00 


12/07/1999 


AntuH 


-FORS1 


2x900 


V 


1.77 


0.54 




1 


13 


26 


23.5 


-42 52 00 


08/04/1999 


AntuH 


-ISAAC 


10x36x10 


K a 


1.15 


0.41 


A 


1 


13 


26 


23.5 


-42 52 00 


07/05/1999 


AntuH 


-ISAAC 


10x36x10 


K a 


1.05 


0.40 


B 


1 


13 


26 


23.5 


-42 52 00 


28/05/1999 


AntuH 


-ISAAC 


10x36x10 


K s 


1.14 


0.51 


C 


1 


13 


26 


23.5 


-42 52 00 


15/04/2000 


AntuH 


-ISAAC 


65x6x10 


K 3 


1.09 


0.35 


D 


1 


13 


26 


23.5 


-42 52 00 


10/05/2000 


AntuH 


-ISAAC 


65x6x10 


K a 


1.22 


0.43 


E 


1 


13 


26 


23.5 


-42 52 00 


07/06/2000 


AntuH 


-ISAAC 


65x6x10 


K a 


1.13 


0.60 


F 


1 


13 


26 


23.5 


-42 52 00 


08/07/2000 


AntuH 


-ISAAC 


65x6x10 


Ks 


1.15 


0.31 


G 


2 


13 


25 


24.0 


-43 10 00 


11/07/1999 


AntuH 


-FORS1 


2x900 


U 


1.28 


0.53 




2 


13 


25 


24.0 


-43 10 00 


11/07/1999 


AntuH 


-FORS1 


2x900 


V 


1.20 


0.44 




2 


13 


25 


24.0 


-43 10 00 


08/04/1999 


Antu+ISAAC 


10x36x10 


K a 


1.05 


0.36 


A 


2 


13 


25 


24.0 


-43 10 00 


07/05/1999 


AntuH 


-ISAAC 


10x36x10 


K a 


1.10 


0.39 


B 


2 


13 


25 


24.0 


-43 10 00 


29/05/1999 


Antu+ISAAC 


10x36x10 


KS 


1.06 


0.44 


C 


2 


13 


25 


24.0 


-43 10 00 


15/04/2000 


Antu-i 


-ISAAC 


65x6x10 


K s 


1.09 


0.40 


D 


2 


13 


25 


24.0 


-43 10 00 


12/05/2000 


Antu- 


-ISAAC 


51x6x10 


K s 


1.22 


0.51 


Eb 


2 


13 


25 


24.0 


-43 10 00 


19/05/2000 


Antu+ISAAC 


20x6x10 


Ks 


1.22 


0.44 


Ec 


2 


13 


25 


24.0 


-43 10 00 


07/06/2000 


AntuH 


-ISAAC 


41x6x10 


Ks 


1.13 


0.57 


Fa 


2 


13 


25 


24.0 


-43 10 00 


09/06/2000 


Antu^ 


-ISAAC 


30x6x10 


Ks 


1.13 


0.40 


Fb 


2 


13 


25 


24.0 


-43 10 00 


08/07/2000 


Antu+ISAAC 


65x6x10 


KS 


1.15 


0.40 


G 


2 


13 


25 


24.0 


-43 09 04 


20/02/2000 


NTT+SOFI 


45x10x6 


K s 


1.04 


0.55 




2 


13 


25 


24.3 


-43 09 58 


27/06/1998 


HST+NIC3 


256 


F222M 









did not use the ESO pipeline reduction product. Good 
sky subtraction in a crowded field like that of a galactic 
halo is particularly important. For that step we used the 
DIMSUM package (Stanford et al. |1995[ ) within IRAFQ. 
In DIMSUM the sky subtraction is made in two passes. 
In the first, a median sky is computed for each image 
from the six frames that are closest in time. The shifts be- 
tween the sky-subtracted frames are then computed and 
all the images stacked together using a rejection algorithm 
to remove cosmic rays. An object mask is computed for 
the coadded image and then shifted back in order to cre- 
ate object masks for the individual frames. In the second 
pass, the sky subtraction is made using the object masks 
to avoid overestimation of the sky level. These masks are 
also used to check that the bright object cores were not 
removed as cosmic rays in the previous pass. After the 
mask-pass sky subtraction, all frames belonging to single 
epoch are registered with imalign and combined with im- 
combine task in IRAF. 



3. The Photometry 

3.1. Data analysis 

The photometric reduction of the combined U-, V- and 
-ftTg-band images was performed using the DAOPHOT II 



programme (Stetson 1987, 1994). First, we located all the 



objects that were > 3er above the background on individ- 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under contract 
with the National Science Foundation 



ual images. More than 50 relatively bright, not saturated, 
isolated, stellar objects were chosen to create the variable 
PSF for each image. ALLSTAR fitting of the PSF to all 
the objects produced the object lists that were matched 
with DAOMATCH and DAOMASTER, where only ob- 
jects with good photometry in at least two frames were 
kept. The final photometric catalogue was obtained with 
ALLFRAME which uses as input information photome- 
try lists from ALLSTAR and fits the PSF to all frames 
(in U, V and X-band) simultaneously. Again, only the 
objects detected in at least 2 images were kept. Using the 
information on the location of stars in all bands simulta- 
neously improved our photometry, which is deeper by ~ 1 
mag with respect to ALLSTAR photometry. Moreover, the 
treatment of close companions, in particular the ones that 
have different colors, is much better with ALLFRAME. 

NGC 5128 is close enough that its globular clus- 
ters appear slightly resolved, in the sense of having a 
larger FWHM and non-stellar PSF (Minniti et al. |1996, 



Rejkuba |200l| ) . Restricting the sharpness and goodness of 
the fit parameters to —0.7 <SHARP< 0.7, we rejected 
most of the galaxies, star clusters and other extended ob- 
jects as well as remaining blemishes and cosmic rays from 
the final photometry list. Also stars with large photomet- 
ric uncertainties in one or more filters (a > 0.5 mag) were 
rejected. Finally the images were visually checked and a 
few remaining extended background objects (e.g. partially 
resolved star forming region in a background spiral galaxy) 
were discarded. With this selection our final U, V photom- 
etry catalogue contains 1581 and 1944 stars in Fields 1 and 
2, respectively, while the V, K catalogue contains 5172 and 
8005 stars and the numbers of stars with good photome- 
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Fig. 1. Combined U- and V- band image of Field 1, taken at the VLT (UT1+FORS1). The field of view is 6.'8 x 6.'8. 
North is at the top and east to the left. 



try in U, V and if-band are 508 and 663 in Fields 1 and 
2, respectively. 



3.2. The photometric calibration 

For the photometric calibration of the optical images, 
standard stars from the catalogue of Landolt (1992) were 
used. We checked the photometric quality of the two 
nights separately and since both were photometric, with 
almost identical zeropoints, we combined the standard 
stars for the two nights. Using a total of 14 stars in 4 dif- 
ferent fields, spanning the color range —1.321 < (U—V) < 
4.162, we derived the following calibration transforma- 



tions: 

Uinst = t/- 24.262(±0.089) + 0.379(±0.068)*X 

-0.042(±0.007) * (U - V) (1) 

Vinst = V - 27.348(±0.042) + 0.213(±0.034) * X (2) 

where X is the mean airmass of the observations, Uinst 
and Vi ns t are instrumental magn itudes and U and V mag- 
nitudes from the Landolt ( 1992 ) catalogue. The one-sigma 
scatter around the mean was 0.031 mag for the U band 
and 0.023 mag for the V band (Fig. ||, left panel). Adding 
the color term (U — V) in the transformations slightly 
reduces the scatter in the U-band to 0.017 mag (Fig. |[ 
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Fig. 2. Combined U- and V- band image of Field 2, taken at the VLT (UT1+FORS1). The field of view is 6.'8 x 6.'8. 
North is at the top and east to the left. 



right panel). However, the calibration equation without 
the color term for the y-band was preferred, since the 
scientific data in that filter go much deeper and U magni- 
tudes for some objects could not be measured accurately 
enough. 

Observations at our reference if s -band epoch (G in 
Table |l|) were taken in photometric conditions. During 
the same night, July 8 200 0, three standard stars from the 
list of Persson et al. ( 1998| ) were observed. Each standard 
star was observed at 5 different positions on the IR-array. 
In this way, a total of 15 independent measurements were 
obtained. However, the number of measurements with dif- 
ferent airmass was only 3, so that we preferred to adopt 



the mean extinction coefficients measured on Paranal for 
K s - band of 0.05 mag/airmass. The derived zeropoint of 
the G-epoch observations is 24.23 ± 0.04. The following 
calibration equation was applied to our data: 



ksMst =K S - 24.23(±0.04) + 0.05 * X 



(3) 



where X is the airmass of observations, k s ^ ns t is the in- 
strumental magnitude and K s is the calibrated magni- 
tude. The photometry of all other X-band epochs were 
measured with respect to the reference epoch. 
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Fig. 4. The completeness and photometric uncertainties 
as a function of magnitude and color. The thin vertical 
dotted lines indicate the value of input U — V colors. The 
full line defines the 90% completeness limit and the dashed 
line is for 50% completeness limit. 



3.3. Completeness and contamination 

We made extensive tests to measure completeness and 
magnitude uncertainties as a function of magnitude and 
radial distance from the center of the galaxy. The com- 
pleteness for the U, V and X s -band photometry has 
been calculated using the ADDSTAR programme within 
DAOPHOT. We made twenty artificial star experiments, 
adding each time ~ 3000 stars to the first frame. The stars 
were added on a regular grid separated by 2 x Rpsf + 1 
pixels in order not to alter crowding (where Rpsf is the 
PSF radius used for fitting the image with the worst see- 
ing) and having magnitudes randomly distributed in the 



100 r' * _ 

« « » e „ . I 
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60 ~r 

40 L o F1 

20 : , }F2 U = 24.5 

6-5 H — i — i — i — F 

» 100 r 71 \ 



o o o 
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« 60 : 



I 40 r o fi 

B 20 r , F2 V = 24.0 
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80 2 
60 2 

40 2 o F1 

20 2 , F2 K=21.0 

^ 1 1 1 — - 



o o _ o 



10 15 

R (arcmin) 

Fig. 6. The completeness as function of radial distance 
from the center of NGC 5128 calculated around 90% com- 
pleteness limit for U- and F-band (for (U— V) = 0.7) and 
around 65% for i^ s -band. 



observed range. The position of the first star in the list 
was chosen randomly, so that over 20 different experi- 
ments the added stars were uniformly distributed over the 
whole field. After the appropriate coordinate shifts were 
applied and magnitudes changed to instrumental system 
taking into account the observed magnitudes and colors, 
the same stars from the first frame were added also to all 
other frames. Their photometry was recomputed in the 
same way as for the original images. The stellar PSF ob- 
tained from the field stars for the respective image was 
used in the simulations. Incompleteness in Field 1, defined 
by a recovery rate of 50% from the artificial-star exper- 
iment, sets in around magnitude 25 in the [/-band (but 
depends strongly on U — V color; see Fig. |) and 22.5 m 
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X s -band (Fig. ||). The corresponding numbers for Field 2 
are 25 for [/-band and 21.3 for if s -band. The difference at 
-ftf s -band between the two fields is due to (1) better seeing 
in Field 1 (FWHM = 2.1 pix vs. 2.7 pix) and (2) higher 
surface brightness in Field 2. 

Dependence of the completeness on radial distance 
from the center of the galaxy («20oo = 13 h 25 m 26!4, 
<520oo = — 43°01'05'.'l) was calculated for the magnitude 
bin around 90% level of completeness in U- and F-band 
(Fig. |^) and around 65% of completeness in i^-band. 
There is no significant spatial variation of completeness 
in our data. 

To assess the accuracy of our photometry, we calcu- 
lated the difference between the input and recovered mag- 
nitudes for each magnitude bin (Fig. |I| ||). Our photome- 
try is reliable down to the incompleteness limit and blend- 
ing is not seriously affecting our data (see the discussion 
in Sec. 4.2). For magnitudes fainter than the 50% com- 
pleteness limit, the measured values are systematically 
brighter, because of the bias towards brighter fluctuations 
and blending due to crowding. The colors of the recovered 
stars with magnitudes fainter than the 50% completeness 
limits are redder than the input colors in the UV CMDs 
(Fig. [!]), due to the larger incompleteness in the U- than in 
the 1^-band. In the VK CMDs the colors of the recovered 
stars range from slightly redder than the input color for 
the very red stars, due to the dominant incompleteness in 
the V-band, to bluer for the faintest and bluest stars, due 
to the dominant incompleteness in the if-band (Fig. 
We did not correct our data for this systematic shift, be- 
cause magnitudes fainter than the 50% completeness limit 
will not be used in further analysis. 

Contamination by foreground Galactic stars and un- 
resolved background galaxies is important because of the 
low Galactic latitude of our fields (b — 19?5) and the very 
deep photometric limits observed. We used the Besancon 
group model of stellar population synthesis o f the Galaxy 
avail able t hrough the WebQ (Robin & Creze 198£, Robin 
et al. 1996 ) to simulate the total number and optical mag- 
nitude and color distribution of Galactic foreground stars 
in our fields. The simulated catalogue has 1827 stars in 
the FORS1 field (6f8 x 6f8) in the magnitude interval 
18 < V < 30. In order to get a realistic estimate for 
the number of stars that would be observed in our fields, 
the correction for completeness is necessary. Therefore, we 
added the stars from the simulated catalogue (with mag- 
nitudes scaled to correspond to instrumental magnitudes) 
to our images and re-measured their magnitudes. In this 
way, realistic photometric uncertainties were applied and 
the number of stars recovered in two fields was corrected 
for completeness. A total of 340 and 350 stars with El- 
and V-band photometry satisfying profile fitting and pho- 
tometric uncertainty selection criteria were measured in 
Field 1 and 2, respectively. Most of the foreground stars 
have < U— V < 3 (Fig. 0). All of the stars in the red part 
of the Field 2 CMD brighter than V ~ 22 are expected to 




Fig. 7. UV (left) and VK (right) color-magnitude dia- 
grams for foreground Galactic stars simulated using the 
Besangon group model. The correction for completeness 
and realistic photometric uncertainties have been applied 
by adding the simulated stars to our Field 2 images and 
re-measuring them. The reddening vector E{B — V) = 0.1 
is plotted in the upper right corner. 



be foreground stars (see Sect. 4.1). In order to adjust for 



the expected number of Galactic stars we normalized the 
models to the observed number of reddest stars in Field 2. 
Thus the total number of foreground stars was increased 
by 31% in both fields. 

The Galaxy model simulation supplies not only the 
colors of the simulated stars, but also their metallicities, 
ages, spectral types an d lum inosity classes. Using all these 
data and the Kurucz ( 1998 ) model atmospheres^, we de- 
rived the foreground contamination in the i^-band. After 
correction for completeness, the expected number of fore- 
ground stars in ISAAC images is 112 and 91 in the VK 
CMD of Field 1 and 2, respectively. All of the foreground 
stars have 1.0 < V — K < 4.5 (Fig. right panel). 

The measured number of compact background galaxies 
on the FORS1 images, taken in similar observing condi- 
tions, is 55 400 in magnitude range V = 20 — 25 mag 
for the selection of sharpness parameter — 1 < sharp < 1 
(Jerjen & Rejkuba 2001). Our tighter selection criteria 
(—0.7 < sharp < 0.7) eliminated most of them. In the 
smaller field of view of ISAAC, the predicted number of 
background galaxies is ~ 330 in the interval of magnitudes 
16 < K s < 23 and ~ 180 between 16 < K s < 22 (Saracco 



et al. 2001). Most of the background galaxies are resolved 



and rejected by sharp and magnitude uncertainty param- 
eter requirements on our photometry. Only few compact 
galaxies might contaminate the sample. 



http:/ /www. obs-besancon.fr/www/modele/modele_ang. html 
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Table 2. Comparison with Fasset & Graham (FG) pho- 
tometry; the numbers of stars in the first column are from 
FG 



# 


Vfg 


(B-V)fg 


(U-B) FG 


V'forsi 


C^FORSl 


1 


19.94 


+0.23 


-0.54 


19.987 


19.498 


4" 


21.07 


+0.12 


-1.26 


20.871 


19.711 


5 b 


21.19 


+0.04 


-1.26 


20.925 


19.717 




21.51 


-0.07 


-0.34 




20.291 


10 


21.65 


+0.04 


-1.16 


20.986 


20.440 


12 


21.72 


+0.21 


-0.20 


21.722 


21.580 


13 


21.75 


+0.11 


-0.84 


21.667 


20.600 


15 


21.86 


+0.09 


-0.61 


21.885 


21.099 


16 


21.94 


-0.05 


-0.92 


21.978 


20.862 


17 


21.95 


+0.34 


-0.82 


21.963 


21.391 


18 d 


22.01 


+0.30 


-0.84 






20 


22.10 


+0.07 


-0.82 


22.150 


21.144 


21 


22.12 


-0.08 


-1.00 


22.127 


20.942 


26 


22.25 


+0.08 


-0.98 


22.323 


21.234 


27 


22.35 


-0.03 


-0.91 


22.331 


21.295 


28 


22.40 


-0.05 


-0.87 


22.412 


21.323 


30 


22.51 


+0.20 


-0.56 


22.631 


22.032 


32 


22.53 


+0.13 


-0.50 


22.579 


22.217 


33 


22.54 


+0.07 


-0.76 


22.623 


21.528 


35 


22.55 


+0.11 


-0.55 


22.503 


21.968 


36 


22.58 


-0.08 


-0.88 


22.606 


21.535 


37 


22.58 


+0.02 


-1.06 


22.777 


21.549 


39 


22.69 


-0.12 


-0.82 


22.706 


21.585 


42 


22.82 


+0.00 


-1.02 


22.856 


21.646 


46 


22.89 


+0.31 


-0.12 


22.919 


23.260 


47 


22.90 


-0.05 


-1.00 


22.934 


21.737 


48 


22.96 


+0.08 


-0.74 


23.051 


22.149 


50 


22.99 


-0.04 


-0.81 


23.134 


22.027 



a fl.GC-8 (Rejkuba |2001|) 
b fl.GC-25 (Rejkuba pOOll ) 
c saturated in FORS1 V image 
d extended in FORS1 images 



3.4. Comparison with published data 

Several recent studies of resolved stellar populations in 
NGC 5128 exist in the literature, all but one made with 
HST in F606W (V) and F814W (I) or F110W (J) and 
F160W (H) photometric bands, which are not very sensi- 
tive to the recent star formation. Our Field 2 is centered 



on the field of Soria et al. ( 1996 ), which was also observed 
in t he nea r-IR (F110W and F160W filters) by Marleau et 
al. ( 2000 ), but the direct comparison is not possible due 
to different photometric bands used. 

The only possible direct comparison is for our Field 1 
photometry, which partially overlaps with the HST pho- 
tometry of Mould et al. ( [2000 ) and the ground based pho- 
tometry of Fasset & Graham (2000). The last authors ob- 
served a wider field than ours in U, B and V Glass fil- 
ters, but due to smaller telescope aperture (2.5m), their 
photometry is much shallower. The HST photometry is 
deeper, but covers a much smaller area. The mean dif- 
ference between V magnitudes of Fasset & Graham and 




Fig. 8. Com paris on of our photometry with that of Fasset 
& Graham (pooop for 26 bright blue stars in common. In 
the upper panned the comparison of V band magnitudes 
is presented and in the lower panel U-band magnitudes 
are compared. The two exte nded objects (#4 and #5 in 
Tab. |) noted by Rejkuba (|200l| ) are plotted with filled 
symbols. The mean offset of -0.13 mag in t/-band is indi- 
cated with the dashed line. 

Mould et al. photometry is —0.13 ± 0.07 mag, in the sense 
of HST photometry having a systematically fainter zero 
point (Fasset & Graham [2000 ). Comparison of 26 stars in 
common between our data and that of Fasset & Graham 
for the brightest blue stars (their Table 3) is presented in 
Fig. ||. The mean difference for all 26 stars is negligible 
for y-band photometry, but it amounts to 0.13 mag in U- 
band. From Fig. ^ a systematic trend with the magnitude 
is apparent in V and U band. Excluding the star #10 (see 
Tab. |), which has HST V magnitude of 21.04 (thus close 
to our value; star R4 in Tabic 1 of Mould et al.), and ex- 
tended sources #4, #5 (see Rejkuba 2001 ; filled triangles 
in Fig. ||) and #18, these trends can be represented by the 
following equations: 



Vtorsi - ^fg = 0.030 x V FO rsi - 0.65 
C/forsi ~ U FG = 0.089 x E/forsi - 2.05 



(4) 
(5) 



with rms equal to 0.055 for V and 0.096 for U. At least 
part of this difference may be due to different filters used 



http : // cfaku5 .harvard .edu / grids, html 



(Bessell 1995). Fasset & Graham further neglected the 
color term in their calibration due to insufficient color 
coverage of their standards, while we found that the color 
term is important for our U -band calibration. 

4. The Color-Magnitude Diagrams 
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Fig. 9. UV color-magnitude diagram for Field 1 (left) and Field 2 (right). The arrow represents the reddening vector 
of E{B — V) = 0.1. Typical photometry uncertainties are plotted on the right side. The solid line indicates the 50% 
completeness. 



4.1. Optical CMDs 

Optical CMDs probe young and intermediate age popu- 
lations. Theoretically, the U and V light is dominated by 
young main sequence stars (Buzzoni 1995). 

Figure ^ shows (U — V) -V color-magnitude diagrams 
for stars in both observed halo fields of NGC 5128. Due to 
better seeing the saturation magnitude of the y-band is 
~ 0.3 mag fainter for Field 2. Most of the stars redder than 
(U — V) ~ belong to our own Galaxy (see Fig. 0). The 
most important characteristic of the UV CMD of Field 1 
is the upper main sequence, visible as the blue plume at 
(U — V) ~ —1 mag. By contrast, there are no such young 
massive stars in Field 2. 



The Besancon Galaxy model (see previous section) was 
used to "clean" the CMDs of foreground star contamina- 
tion. In Fig. |l^ we show UV CMDs after the subtraction of 
foreground stars. Overplotted are isochrones from Bertelli 



et al. Q1994D for Z = 0.004 and for log(age)=7.0, 7.3, 7.5 
for Field 1 and log(age)=7.0, 7.5 and 7.8 for Field 2. The 
isochrones were shifted to the distance of NGC 5128, as- 
suming a distance modulus of (m — M)v = 27 .8 an d red- 
dening of E(B-V) =0.11 mag (Schlegel et al. |l998[) . The 
difference between the two diagrams is striking: the blue 
main sequence containing stars as young as ~ 10 Myr that 
is present in Field 1 is completely absent from the CMD of 
Field 2. Well separated from the main sequence in Field 1 
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Fig. 10. CMDs with foreground stars statistically subtracted. Overplotted are isochrones from Bertelli et al. (1994) 
for Z = 0.004 with log(age) indicated with the number on the right side. The size of the reddening vector corresponds 
to E(B -V) = 0.1 mag. 



is the sequence of blue core-helium burning (BHeB) stars 
at (U — V) ~ 0.5. The width of the gap and the tightness 
of the main sequence indicate low differential extinction 
in the field. In Field 2, the stars with 0.5 < (U - V) < 2.8 
and brighter than V ^ 22.5 have no corresponding main 
sequence stars and thus cannot be young stars migrating 
from blue to red during their core-He burning phase, as 
is the case for most of the objects in Field 1 with colors 
0.5 < (U — V) < 2.5. These stars in Field 2 are most prob- 
ably the remaining foreground contamination, indicating 
that the foreground contamination may affect the num- 
bers of blue and red HeB stars in Field 1 (see next para- 
graph). Only a few stars are lying along the isochrone of 



log(age) = 7.5 (Fig. |l0|), while most of them are consistent 
with much older ages. We conclude that there are no stars 
younger than ~ 40 Myr in Field 2. 

Note that isochrones for metallicities higher than Z = 
0.004 extend on the red supergiant edge to redder val- 
ues of U — V than the reddest stars in Field 1 and thus 
do not fit well our observations. The ratio of blue to red 
supergia nts s trongly depends on metallicity (Langer & 
Maeder |l995i Maeder & Meynet fcOOlj) and in principle 
could be used to constrain the metallicity of the youngest 
population in NGC 5128. Counting the number of blue 
and red supergiants for stars more massive than ~ 12 M Q , 
we find their ratio to be B/R < 0.7 — 0.8 (although with 
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high uncertainty due to the possible foreground contami- 
nation), in good agreement with the observed B/R value 
in SMC clus ter NGC 300 (see discussion by Langer & 
Maeder |1995| ). The metallicity of Z = 0.004 (correspond- 
ing to [Fe/H]= —0.7 dex) is appropriate for the SMC. 
However, since the B/R value depends also on other pa- 
rameters such as stellar mass, rotation and degree of over- 
shooting (Maeder & Meynet 2001), a more detailed com- 
parison with models and the determination of the metal- 
licity of this youngest stellar population in NGC 5128 is 
warranted (Rejkuba et al., in preparation). 

The low metallicity implied by the fit of the isochrones 
on UV CMDs probably reflects the metallicity of the gas 
left in the halo of NGC 5128 by the accreted satellite. 
Atomic Hi (Schiminovich et al. |1994 ) and molecular CO 



gas (Charmandaris et al. 2000) present in Field 1 are 



slightly offset from the position of the diffuse stellar shell, 
the obvious remnant from the accreted galaxy. 

4.2. Optical-Near IR CMDs 

Optical-near IR CMDs probe old and intermediate-age 
stellar populations. Theoretically, more than two thirds of 
the light in A-band is dominated by cool stars on the red 
giant branch (RGB) and as ympto tic giant branch (AGB), 
and by red dwarfs (Buzzoni 1995 ). The red dwarfs are too 
faint to be detected at the distance of NGC 5128, and thus 
our VK CMDs are entirely dominated by RGB and AGB 
stars (Fig. |ll|). 

In Fig. H we show VK CMDs of Fields 1 and 2 after 
the statistical subtraction of foreground stars. Overlaid 
are fiducial RGB sequences of Galactic globular clusters 
(from left to right: M15 , 47 Tuc, NGC 6553 and NGC 
6528; Ferraro et al. |2000| ) spanning a large range of metal- 
licities (-2.17 <[Fe/H]< -0.23 dex). As before, we used 
a distance modulus of 27.8 and reddening corresponding 
to E{B - V) = 0.1 (E(y - A") = 0.274 and A K = 0.0347; 
Ricke & Lebofsky 1985) to adjust the magnitudes and 
colors of RGB fiducials to those of NGC 5128. Obviously, 
most of the stars in Fig. [l^ belong to the RGB. The right 
edge of the RGB is quite sharp, with most of the stars be- 
ing more metal-poor than 47 Tuc ([Fe/H]= —0.71 dex) 
and none appearing to be as metal-rich as NGC 6553 
([Fe/H]= —0.29 dex). However, the latter is due to in- 
completeness in V-band photometry. 

The spread in color of the RGB is larger than the 
photometric uncertainties (Fig. |l3|), indicating the pres- 
ence of spread in metallicity and/or age. The most metal- 
poor stars have metallicities of —2 dex if their ages cor- 
respond to those of Galactic globular clusters. The popu- 
lation we probe with VK photometry is more metal-poor 
than —0.7 dex. Our V-band images are not d eep e nough 
to detect more metal rich giants. Walsh et al. ( 1999 ) mea- 
sured the mean oxygen abundance of five planetary nebu- 
lae in NGC 5128 to be [0/H]= -0.5 ±0.3 dex, consistent 
with the presence of the large population of stars with 
metallicities below solar, as we observe in the VK CMDs. 
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Fig. 13. The histogram of the observed color distribution 
along the RGB (lower panels) is plotted accounting for 
the actual slope of the RGB (thick line in upper panels) 
in Field 1 (left) and 2 (right). The lcr spread in color of 
the RGB is indicated in the upper left corner. Overplotted 
with the dashed line is a Gaussian with the a — 0.3, typical 
of the distribution from the photometry errors. 

There are 1830 stars in Field 1 and 1197 in Field 2 
which have good A-band photometry [uk < 0.5, —0.7 < 
sharp < 0.7 and \ < 2.0) above the respective A-band 
completeness limits, but which have not been detected in 
optical bands. They are uniformly distributed over the 
whole ISAAC images. The large number of very red stars 
with good photometry in K s and no counterpart in V- 
band suggests that stars more metal- rich than [Fe/H] = 
—0.7 dex are present, as expected for a luminous giant 
elliptical galaxy. This is in good agreement with the results 
of Harris et al. Q1999Q and Harris & Harris ( |2000D . 



5. The Color-Color Diagram 

Fig. [u] shows the color-color diagram for the objects de- 
tected in U-, V- and A-band in Field 1 (open squares) and 
Field 2 (filled triangles) that had ALLFRAME photome- 
try uncertainties < 0.5 mag and magnitudes brighter than 
the 50% completeness limits in all bands. The error-bars 
plotted are the ones given by ALLFRAME. A total of 99 
objects in Field 1 and 90 in Field 2 satisfy these criteria. 
Most of them are located along the stellar sequence that 
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Fig. 11. VK color-magnitude diagram for Field 1 (left) and Field 2 (right). The dashed line identifies the 50% 
completeness level. 



crosses the color-color diagram in diagonal. Matching the 
UV and VK catalogues after the statistical subtraction of 
the foreground contamination, the stellar sequence disap- 
pears and only 54 and 44 objects with good photometry 
in all three filters are left in Field 1 and 2, respectively. 

We use the color-color diagram presented in Fig. [li] 
to measure the reddening and the contamination by back- 
ground galaxies. 

5.1. Reddening 

The objects crossing the color-color diagram in diagonal 
from top right to bottom left belong to the Milky Way. The 
quality of the fit of the stellar sequence with isochrones de- 
pends on the age of the isochrone with respect to stars. 
Because the foreground Milky Way stars span a range of 
ages, isochrones of different ages fit different parts of the 
stellar sequence. As an example, in order to guide the eye 
we overplotted the stellar loci of solar metallicity stars 
and log(age)=7.0 yr (Bertelli et al. |l994| ). The tight stel- 
lar sequence in the U — V vs. V — K color-color diagram 
shows that the foreground reddening in the two fields is the 
same. Its magnitude was measured by fitting the theoret- 



ical isochrones and it amounts to E{B — V) = 0.15 ± 0.05 
ma g. Thi s is in excellent agreement with the Schlegel et 
al. ( 1998 ) extinction maps a nd the reddening determined 
by Fasset & Graham ( 200C| ), who used UBV color-color 
diagrams to measure a value of E(B — V) = 0.14 ± 0.02 
mag. 

5.2. Star-galaxy separation 

Using the color-color plots it is possible to distinguish the 
foreground dwarf stars from the most compact unresolved 
galaxies. Galaxies are located in the lower right part of the 
diagram in Fig. |lj. The resolved galaxies are not plotted 
here since they got rejected by our selection of shape pa- 
rameters in FIND and ALLFRAME. Note that the group 
of objects with —0.5 < (V — K) < 2 are the young stars 
detected in Field 1. 



6. Jet aligned Star Formation in the Halo 

The chains of bright blue compact obj ects a nd optical fil- 
aments, first noticed by Blanco et al. ( 1975 ), were recog- 
nized as young blue supergiants and OB associations by 
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Fig. 12. VK CMDs with foreground stars statistically subtracted. Overplotted are fiducial RGB sequences of Galactic 
globular clusters (from left to right): M15, 47 Tuc, NGC 6553 and NGC 6528, from Ferraro et al. (§00(3). 



Graham & Price (fL98l|) and Graham ( |l998|) . Our Field 1, 
situated ~ 14 kpc away from the center of the galaxy and 
coinciding with the NE radio lobe, contains part of this 
"string of blue knots and filaments." 

In Fig. [l5| we show the spatial distribution of the bluest 
(left) and the reddest (right) stars in both our fields. The 
red stars, belonging to the foreground Galactic popula- 
tion, are uniformly distributed in both fields. Field 2 has 
almost no very young blue stars, while the majority of 
the young population in Field 1 is aligned with the jet 
direction coming from the nucleus of the galaxy. The ver- 
tical structure to the right of the jet is located at the left- 
most edge of a large Hi cloud found by Schiminovich et 
al. ( |1994j ) . In the northern part of this field, a CO molec- 
ular cloud was recently discovered by Charmandaris et 
al. ( |2000D . 

Mould et al. (|2000D, Graham (p99Sh and Fasset & 



Graham (2000) report on the presence of the collimated 
star formation in the halo of NGC 5128. They favor the 
scenario in which a past i nteraction between the radio 
jet (Morganti et al . |l99S| ) and the Hi cloud complex 



(Schiminovich et al. 1994[) is responsible for the star for- 
mation. 



Finding active star formation far out in the halo of an 
elliptical galaxy is unusual. Understanding the triggering 
mechanism and the origin of the star formation, as well as 



the nature of the underlying ionized gas (Graham 1998 



Morganti et al. 1991, 1992), is important for explaining the 



filamentary emission and jet-aligned structures observed 
in other galaxies (e.g., in NGC 1275: Lazareff et al 



more distant radio sources: Rees 1989, Best et al. 1996, De 



Breuck et al. 1999) 



1989 



Sabra et al. |200C| ; in M87: Gavaz zi et a l. |2000| ; an d in ot her 



The alignment of the luminous blue stars with the ra- 
dio axis of NGC 5128 (Fig. [l5]) may shed some light on 
the debated issue of the observed alignment of the radio 



and optical structures in high-z radio galaxies (Rees 1989 
Best et al. 1996). One of the proposed alignment mod- 



els is based on the assumption that the action of the jets 
tunneling through a galactic medium affects the gas via 
shock heating, ther eby in ducin g vig orous star f ormation 
(Begelman & Ci ofii |l989| , Rees |l989| , Daly |l990|) . In par- 
ticular, as Rees (198E) discussed in detail, the stars form 
from the cool-phase gas (with T < 10 4 K) in clouds mas- 
sive or dense enough to be Jeans-unstable. 
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Fig. 14. Color-color diagram for the objects matched in 
U, V and K s frames that had ALLFRAME photometry 
uncertainties < 0.5 mag and magnitudes brighter than the 
50% completeness limits in all bands. Field 1 objects are 
plotted with open squares and Field 2 with filled tria ngles . 
Overplotted is the Padova isochrone (Bertelli et al. 1994 ) 
for the solar metallicity and log(age)=7.0 yr that indicates 
the locus of the foreground stellar sequence. Reddening 
vector has dimensions corresponding to E(B — V) =0.1. 



This mechanism may work well for massive high-z ra- 
dio galaxies, in which most of the baryonic mass would 
still be gaseous. In the present day giant ellipticals, on 
the other hand, little gas is usually left to support strong- 
star formation. However, probably due to previous accre- 
tion of a gas rich galaxy, NGC 5128 has substantial cool, 
dense clouds of mole cular CO and atomic Hi gas in its halo 
(Schiminovich et al. 1994 , Charmandaris et al. 2000|) . We 
note that these clouds could form the clumpy ISM/IGM 
necessary for the 'bursting bubble' model proposed by 
Morganti et al. ( |1999| ). 

We have shown clear evidence for blue stars as young 
as 10 x 10 6 yr (Fig. |hJ). We also see underlying emission, 
particularly in our £/-band images, which probably comes 
from the [Oil] 3727 emission line. High-resolution spectra 
with good signal-to-noise are needed to decide whether 
this emission is due to gas being photoionized by the newly 
formed stars or being photoionized by particles coming 
directly from the nucleus. 

Other low-redshift galaxies also show filamentary ion- 
ized gas struct ures in their halos. In the case of M87 
(Gavazzi et al. 



2000[ ), the filament of ionized gas in the 



NE part of the halo coincides with the Eastern radio-lobe, 
in much the same way as in the much closer NGC 5128. It 
would be interesting to search for young stars associated 
with that filament too. 



Fig. 15. The spatial distribution of the bluest (U — V) < 
—0.75 (left panels) and the reddest (U — V) > 2 (right 
panels) stars in Field 1 (top panels) and 2 (bottom panels) 
for stars with magnitudes above the 50% completeness 
limit (V < 25 and U < 25). The north is up and east is 
to the left. 



7. Intermediate Age Population 

In contrast with the conclusions of HST studies (Harris et 
al. 1999, Harris & Harris 2000) of stellar populations in 
NGC 5128, we detect not only old population II stars, but 
also a significant number of stars with magnitudes brighter 
than the tip of RGB {Mj} v = -0.64(±0.12)[M/H] - 



6.93(±0.14); Ferraro et al. |2000J) . 

W e took bolometric corrections from Bessell & Wood 
( 1984 ) and the e mpiri cal fit of (V—K) vs. T c g from Bessell, 
Castelli & Plez fll998| ) to transform our K-(V-K) CMDs 
to the theoretical plane (Fig. O). Overplotted on th e H-R 



diagrams are the Padova tracks from Girardi et al. ( 2000 ) 
for the first-ascent giant branch stars with masses M=0.6, 
0.8, 1.2 and 1.6 M and metallicities Z = 0.004 (full lines) 
and Z = 0.008 (dashed lines). The sharp cut-off on the 
right side of the H-R diagrams is due to incompleteness in 
y-band. 

In discussing the H-R diagram, we should con- 
sider blending. According to theoretical predictions 
(Renzini 1998| ), at the fiducial galactocentric distances of 
9, 12.8 and 17.9 kpc, the number of blends consisting of 2 
stars belonging to the tip of RGB is ~ 1300, 300 and 55 
stars (per 2.2' x 2.2'; i.e. field of view of ISAAC), respec- 
tively. In this calculation, the B-band surface brigh tness 
measurements from Mathieu, Dejonghe & Hui ( 1996 ) were 
used. Obviously, the inner most regions of the galaxy, at 
the distance < 10 kpc from the nucleus are too crowded 
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to give accurate photometry for RGB stars and most of 
the stars above the tip of RGB. 

We attempted to use the NIC3 HST images overlap- 
ping with our Field 2 (galactocentric distance R gc ~ 9 
kpc) data to improve the resolution. However, we found 
out that many faint stars visible on ISAAC images are 
completely within the noise of NIC3 data. The resolution 
of our ISAAC data in the best seeing is as good as that 
of NICMOS images, with the advantage of having much 
better S/N and a more stable PSF. The confirmation of 
the number of AGB stars in Field 2 can be obtained only 
through the analysis of the LPVs (Rejkuba et al. in prepa- 
ration) . 

In Field 1 {R gc ~ 14 kpc), our outer shell field, crowd- 
ing is not that severe. The expected number of two-star 
blends at the RGB tip in this field ranges from 300-50, 
as the surface brightness drops across the field. The total 
number of such blends is therefore less than ~ 200 in the 
most pessimistic calculation. On the other hand, the num- 
ber of stars above the tip of RGB (Afboi ~ —3.7; Girardi 
et al. 200C| ) is 768. Subtracting the number of possible 
blends and allowing for a few of the brightest and bluest 



stars (~ 30 stars with Mr < 9) to be the remaining fore- 
ground contamination, there are still more than 500 stars 
whose position in the H-R diagram and CMD is consistent 
with an intermediate-age AGB population. 

In the inner halo field, Field 2, the number of stars 
above the tip of RGB is 2844. After subtracting the num- 
ber of expected blends 1500), the number of AGB 
stars is 2.5 times larger than in Field 1. This confirms the 
presence of gradients in the intermediate-age population 
within the halo of NGC 5128, as suggested by Marleau et 
al. p00l ). 

The intermediate-age AGB population could have eas- 
ily been missed in V and /-band HST studies due to small 
field of view and small (< 0.5 mag) optical magnitude dif- 
ference between the tip of the RGB and the tip of AGB. 
The AGB stars are up to ~ 2 magnitudes brighter in Mk 
and Mboi than RGB stars and thus are easily detectable in 
near-IR. Thanks to this, Marleau et al. (200C) could detect 
some AGB stars in the much smaller NICMOS field. 

The brightest stars in M32, in the Galactic bulge 
and the bulge of M31 have simila r brightnesses, reach- 
ing Afboi = -5.5 (Freedman [l992|, Elston & Silva |1992, 
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Fig. 17. if-band (left) and bolometric magnitude lumi- 
nosity functions (right) for stars redder than (V — K s ) > 2 
in Fields 1 and 2. 



Frogel & Whitford |1987| , Rich & Mo uld |1991| ). Due to this 
similarity, Davidge & van den Bergh ( 2001 ) suggested that 
the tip of AGB could be used as standard candle for deter- 
mination of distances. In Fig. [l^ we present the J^-band 
and bolometric magnitude luminosity functions. The tip 
of the AGB is observed at bolometric magnitude of —5 
in both fields (in spite of the crowding in Field 2), con- 
sistent with the adopted distance modulus for NGC 5128. 
Most of the stars brighter than this magnitude (as well as 
stars brighter than Mk = — 9; see left panels in Fig. |l7|) 
probably belong to our Galaxy. 



8. The Shell 

In our deep V-band image of Field 1, the prominent NE 
shell is easily recognizable as the region with higher sur- 
face brightness, and also with higher density of stellar ob- 
jects. All the stars that are spatially coincident with the 
shell on the V-band image are shown in Fig. |l8|. In the left 
two panels we plot the UV CMD of all stars coincident 
with the position of the diffuse shell (middle panel), while 
the corresponding VK CMD is on the right. The much 
smaller field of view of ISAAC is centered on the shell so 
that the shell stars cover most of the image. Note that 
projected on the shell is the stream of blue, young objects 
that are aligned with the jet direction. Excluding the stars 
that are coincident with the star forming area (lower pan- 
els), the stellar population belonging to the shell is very 
similar to the one observed in the halo of NGC 5128, i.e., 
our Field 2 population (Fig. Q |ll|). As in the rest of the 
galaxy there is a significant spread in metallicity. 



The presence of stellar shells or ripples around early- 
type galaxies is usually explained as evidence of the accre- 
tion and later merging of a smaller companion galaxy. The 
simulations of the stellar component in the merging galax- 
ies have shown that the shells are created as a result of 
phase-wrapping or of spatial-wrapping of the tidal debris 
(Quinn 1984; Hernquist & Quinn 1989). The companion 
galaxies may contain different stellar populations with re- 
spect to the giant elliptical (e.g. most dwarf galaxies have 
much lower metallicity than giants Mateo 1998). 

Here for the first time we obtained the magnitudes 
and colors of stars belonging to a shell. They are surpris- 
ingly similar to the rest of the halo stars in NGC 5128. 
However, due to incompleteness of our V-band photome- 
try, we cannot probe the most metal rich part of the halo 
of NGC 5128. Very high resolution and deeper observa- 
tions of the shell stars are possible with the Advanced 
Camera for Surveys at Hubble Space Telescope. In order 
to be able to better separate the stars that genuinely be- 
long to shells from the "normal" halo stars in NGC 5128 
and confirm the present result a more quantitative com- 
parison between the stellar populations in our two fields 
is necessary (Rejkuba et al., in preparation). 

9. Summary 

Using the VLT with FORS1 and ISAAC, we have resolved 
stars in the halo and in the diffuse north-eastern shell of 
the closest giant elliptical galaxy NGC 5128. With the 
U , V and iVrband photometry we probe different stel- 
lar populations. The U filter, in particular, is sensitive to 
the youngest stars. It has revealed the striking difference 
in the star formation history between the the inner halo 
field (Field 2) and outer halo field (Field 1), the latter 
coinciding with the diffuse shell that is presumably the 
signature of a recent merger. 

In Field 1, stars as young as 10 Myr are present, while 
there are no stars younger than at least 40 Myr in Field 2. 
Thanks to the high quality of our UV photometry, we 
also detect the gap between the main sequence and the 
blue core-helium burning supergiants. The foreground and 
background contamination has been corrected using the 
Besangon Galaxy model and our own color-color diagram. 

The recent star formation in Field 1 is approximately 
aligned with the direction of the jet coming from the center 
of AGN. However, there is also a chain of young blue stars 
offset from the jet direction that is aligned with the edge 
of the molecular clouds. The presence of the atomic and 
molecular gas (Schiminovich et al. 1994, Charmandaris 



et al. 2000) and their association with the diffuse stellar 
shell is compatible with the dynamical scenario of phase- 
wrapping following the merger of a smaller gas-rich galaxy 
with NGC 5128 (Quinn |1984[ ). In this scenario, the inter- 
action of the material coming from the AGN in the center 
of the galaxy with the material left from the merger stim- 
ulated the star formation in the halo, ~ 15 kpc away from 
the galactic centre. The metallicity of the newly formed 
stars gives an upper limit on the abundance of the gas 
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Fig. 18. The color-magnitude diagrams for the stars located at the position of the diffuse NE shell. 



of the accreted satellite with value of Z = 0.004, a value 
typical of the SMC. 

While the UV CMDs probed the youngest stellar pop- 
ulations, the combination of optical and near-IR filters is 
most sensitive to old and intermediate-age populations. In 
our VK CMDs we observe a very wide and extended giant 
branch. The width indicates the presence of Population II 
stars with metallicities ranging from —2 to —0.7 dex. 
Unfortunately, our V^-band images are not deep enough 
to detect more metal rich stars, but the large number of 
stars detected on if-band images without a counterpart in 
V-band suggests an even redder and more metal-rich RGB 
population, in agreement with the previous HST studies 
(Soria et al. |1996| , Harris et al. |l999[ Harris & Harris |2000D . 

In contrast with the previous two studies, we find a 
very extended giant branch up to Mb ] = —5, which re- 
veals the presence of the intermediate-age AGB popula- 
tion. Our extensive crowding experiments and theoreti- 
cal predictions of the amount of blending (Rcnzini 1998 ) 
demonstrate that blending due to crowding is not signif- 
icant enough to mimic this population in the outer field 
(Field 1). In the inner field, the extent and the properties 
of the AGB population will be assessed through the study 
of the long period variables that have been detected in 
both fields (Rejkuba et al. in preparation). 
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